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Abstract 


The USCLIVAR working group on drought recently initiated a series of global climate model 
simulations forced with idealized SST anomaly patterns, designed to address a number of 
uncertainties regarding the impact of SST forcing and the role of land-atmosphere feedbacks on 
regional drought. Specific questions that the runs are designed to address include: What are the 
mechanisms that maintain drought across the seasonal cycle and from one year to the next? 

What is the role of the leading patterns of SST variability, and what are the physical mechanisms 
linking the remote SST forcing to regional drought, including the role of land-atmosphere 
coupling? The runs were carried out with five different atmospheric general circulation models 
(AGCMs), and one coupled atmosphere-ocean model in which the model was continuously 
nudged to the imposed SST forcing. This paper provides an overview of the experiments and 
some initial results focusing on the responses to the leading patterns of annual mean SST 
variability consisting of a Pacific El Nino/Southem Oscillation (ENSO)-like pattern, a pattern 
that resembles the Atlantic Multi-decadal Oscillation (AMO), and a global trend pattern. 

One of the key findings is that all the AGCMs produce broadly similar (though different in 
detail) precipitation responses to the Pacific forcing pattern, with a cold Pacific leading to 
reduced precipitation and a warm Pacific leading to enhanced precipitation over most of the 
United States. While the response to the Atlantic pattern is less robust, there is general 
agreement among the models that the largest precipitation response over the U.S. tends to occur 
when the two oceans have anomalies of opposite sign. That is, a cold Pacific and warm Atlantic 
tend to produce the largest precipitation reductions, whereas a warm Pacific and cold Atlantic 
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tend to produce the greatest precipitation enhancements. Further analysis of the response over 
the U.S. to the Pacific forcing highlights a number of noteworthy and to some extent unexpected 
results. These include a seasonal dependence of the precipitation response that is characterized 
by signal-to-noise ratios that peak in spring, and surface temperature signal-to-noise ratios that 
are both lower and show less agreement among the models than those found for the precipitation 
response. Another interesting result concerns what appears to be a substantially different 
character in the surface temperature response over the U.S. to the Pacific forcing by the only 
model examined here that was developed for use in numerical weather prediction. The response 
to the positive SST trend forcing pattern is an overall surface wanning over the world’s land 
areas with substantial regional variations that are in part reproduced in runs forced with a 
globally uniform SST trend forcing. The precipitation response to the trend forcing is weak in 
all the models. 

It is hoped that these early results, as well as those reported in the other contributions to this 
special issue on drought will serve to stimulate further analysis of these simulations, as well as 
suggest new research on the physical mechanisms contributing to hydroclimatic variability and 
change throughout the world. 
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1. Introduction 


In recognition of the profound societal impact of drought in many regions of the world and the 
emerging capabilities in simulating drought with global climate models, the USCLIVAR program 
initiated a drought working group in 2006 to “facilitate progress on the understanding and 
prediction of long-term (multi-year) drought over North America and other drought-prone regions 
of the world, including an assessment of the impact of global change on drought processes” 
(Gutzler and Schubert 2007)”. The specific tasks of the working group were to 1) propose a 
working definition of drought and related model predictands of drought, 2) coordinate evaluations 
of existing relevant model simulations, 3) suggest new experiments (coupled and uncoupled) 
designed to address outstanding uncertainties in the nature of drought, 4) coordinate and encourage 
the analysis of observational data sets to reveal antecedent linkages of multi-year droughts and 5) 
organize a community workshop to present and discuss the results. 

This paper provides an overview and some results of task (3) of the working group involving the 
design, coordination, implementation, and initial evaluation of a new set of model simulations that 
address the roles of sea surface temperature forcing and land-atmosphere feedbacks in the 
development and maintenance of drought. This work extends and builds upon recent modeling 
studies (e.g., Hoerling and Kumar 2003, Schubert et al. 2004a and b, Seager et al. 2005, Wang et al. 
2008; Wang et al. 2009) as well as numerous observationally-based studies (e.g., Trenberth and 
Guillemot 1996; Mo et al. 1997, Ting and Wang 1997; Nigam et al. 1999; Koster et al. 2003; Ruiz- 
Barradas and Nigam, 2004; McCabe et al. 2004; Wang et al. 2006) that have provided substantial 
insights into the nature of drought and the important role of both the oceans and land-atmosphere 
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interactions. In particular, this work addresses remaining uncertainties regarding the nature of the 
physical mechanisms linking remote SST forcing to regional drought, the relative contributions of 
the different ocean basins and different time scales of SST variability, and the strength of land- 
atmosphere feedbacks, and thereby starts to frame fundamental questions about the predictability of 
long-term drought. Specific questions addressed by the working group include: What are 
mechanisms that maintain drought across the seasonal cycle and from one year to the next? What is 
the role of the different ocean basins, including the impact of El Nino/Southem Oscillation (ENSO), 
the Pacific Decadal Oscillation (PDO), the Atlantic Multi-decadal Oscillation (AMO), and wanning 
trends in the global oceans? What is the role of the land? To what extent can droughts develop 
independently of oceanic variability due to year-to-year memory that may be inherent to the land? 

In order to address these questions, the working group proposed that a number of mechanistic 
experiments be performed that are designed to address some of the key issues outlined above using 
several different global climate models. A key objective was to be able to assess unambiguously 
the model dependence of the results. In order to accomplish that, it was proposed that each model 
be forced with the same set of idealized SST forcing anomalies. In addition, it was proposed that a 
control run be produced in which each model was forced with the same climatological SSTs. In 
order to allow an assessment of land-atmosphere feedbacks, an additional set of runs was proposed 
in which the soil moisture was fixed using a common approach that could be easily implemented in 
each of the models. The main SST forcing patterns, the experiments, and the models are described 
in Section 2. Section 3 presents some basic comparisons of the model responses to the leading 
patterns of SST variability, with a focus on the United States. The summary and conclusions are 
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given in Section 4. Information on auxiliary experiments and data availability are given in the 
Appendix. 

2. The SST Forcing Patterns, Experiments, and Models 

a) SST forcing patterns 

The basic SST data used in this study are the 1901-2004 monthly SST data produced by Rayner et 
al. 2002). The leading patterns of SST variability are isolated using rotated empirical orthogonal 
functions (REOFs), where VARIMAX rotation (e.g., Richman 1986) is used to help separate the 
leading patterns of Pacific and Atlantic SST variability. 

The REOFs are computed from global gridded values of annual mean SST for the period 1901- 
2004. The use of annual means is meant to address the basic question of the nature of the forcing of 
regional hydroclimates on interannual time scales. This of course, does not distinguish between, 
for example, ENSO and Pacific decadal variability (e.g., Barlow et al. 2001) - the responses to 
leading SST anomaly patterns that occur on these different time scales is addressed with other 
supplementary REOFs described in the Appendix. 

Figure 1 shows the three leading REOFs and associated principal components (PCs) of the annual 
mean SST. The first is a global trend pattern, with wanning occuning over most of the global 
ocean. In fact, this spatial pattern can essentially be reproduced by plotting the slopes of the linear 
trends fit to the 1901-2004 annual mean SST at each grid point (not shown). The associated PC 
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shows a somewhat more complicated (non-linear) long-term time evolution, showing an almost 
step-like increase at about 1940, and a clear trend occurring only after the middle 1960s 1 . The 
second REOF is a pan-Pacific ENSO-like pattern that includes a weak Indian Ocean component. 
The associated PC shows that this pattern varies on both interannual (ENSO) and decadal time 
scales, with the latter including the well known shift that occurred in the middle 1970s (e.g., 
Trenberth and Hurrell 1994). The third REOF is confined for the most part to the North Atlantic 
Ocean and resembles the Atlantic Multi-Decadal Oscillation pattern (Enfield et al. 2001). The 
associated PC shows that while this pattern has clear decadal variability it also exhibits 
considerable interannual variations. These three leading rotated EOFs will in the following be 
referred to as the Trend , Pacific and Atlantic forcing patterns. 

b) Proposed Experiments 

The proposed experiments consist of 50-year 2 simulations in which the model is forced with one or 
more of the idealized SST anomaly patterns (the Trend , Pacific and Atlantic forcing patterns 
described above 3 ). The full forcing patterns are produced by adding scaled versions of the REOFs 
to the long-term monthly-varying SST climatology (defined for the period 1901-2004). The scaling 
factor for the Pacific and Atlantic consists of either plus or minus 2 standard deviations of the 
associated PCs. This rather large amplitude was chosen to help isolate what in some cases may be 
rather subtle SST-drought linkages from relatively short model integrations. In the case of the 
Trend , the pattern is scaled by plus or minus 1 standard deviation, which effectively forces the 

1 We note that there is considerable sensitivity of the trend pattern to the time period of interest and the datasets 
used particularly regarding the contribution from the Pacific (e.g., Vecchi et al. 2008: Chen et al. 2008) 

2 The NOAA Global Forecast System Model (GFS) experiments were somewhat shorter (35 years). 
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model with the 40-year averages of the trend anomalies at the beginning (- 1 standard deviation) 
and end (+ 1 standard deviation) of the 1901-2004 time period. An additional experiment was 
proposed to force the models with a globally-uniform SST wanning distribution equal to the global 
mean of the positive trend pattern (0.1 6°C). 

It is important to note that the anomaly patterns are fixed in time and therefore do not have an 
annual cycle. The absence of a seasonal cycle in the forcing arguably diminishes the importance of 
winter SST anomalies relative to summer SST anomalies, especially for tropical Pacific variability, 
because observations indicate that winter anomalies tend to be much larger in magnitude compared 
to summer anomalies. There is however an annual cycle in the full SST forcing fields as a result of 
the annual cycle in the climatological SST on which the prescribed anomaly patterns are 
superimposed. The SST forcing is repeated with no interannual variability for each year of each 
experiment, but the models still generate interannual variability due to unforced "weather noise" 
associated with the internal dynamics of the models. We will assess the magnitude of the forced 
response to prescribed SST anomalies relative to the magnitude of the unforced interannual 
variability. 

The nine baseline experiments of the project consist of the runs in which the models are forced with 
all eight combinations of the Pacific and Atlantic patterns shown in Fig. 1, as well as the control run 
forced only with the monthly varying climatological SST. These are summarized in Table 1 
according to the different combinations of patterns and phases of the forcing. For example, PwAc 
indicates that a model is forced with the warm phase of the Pacific and the cold phase of the 

3 Masks were applied to REOFs two and three to zero-out any small values that, for the Pacific pattern, fall outside 
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Atlantic patterns. In addition to those shown in Table 1, mils were proposed in which the models 
are forced with either the positive or negative Trend pattern, both alone, or superimposed on 
selected combinations of the Pacific and Atlantic patterns. 

A number of other auxiliary experiments were proposed to isolate further various mechanisms and 
time scales of variability. Some isolate the role of the tropics, while others attempt to separate the 
contributions from ENSO and lower frequency Pacific variability. Another set of experiments were 
formulated to assess the impact of land-atmosphere feedbacks. These additional experiments are 
described in the Appendix. Another related and important set of experiments consist of AMIP-style 
simulations (Gates et al. 1999). These are simulations (typically several decades long) in which the 
models are forced by the historical record of observed SSTs. While the focus of the runs described 
here is on understanding mechanisms and model sensitivity to idealized SST forcing, the AMIP 
runs are important in that they facilitate model validation by allowing more direct comparisons of 
the results with observations. 


Table 1: The different combinations of the Pacific and Atlantic SST anomaly patterns used to 
force the GCMs. Here w refers to the warm phase of the pattern (with a 2 standard deviation 
weight) and c refers to the cold phase (with a 2 standard deviation weight). Also, n denotes 
neutral indicating that the pattern has zero weight. In particular, the PnAn experiment denotes 
the control mil forced with the annually-varying climatological SST. 



Warm Atlantic 

Neutral Atlantic 

Cold Atlantic 

Warm Pacific 

PwAw 

PwAn 

PwAc 

Neutral Pacific 

PnAw 

PnAn 

PnAc 

Cold Pacific 

PcAw 

Pc An 

PcAc 


c) Contributing groups and models 


the Pacific and Indian Oceans, and that, for the Atlantic pattern, fall outside the Atlantic Ocean. 
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A number of groups have participated in this project by contributing model runs. While only a 


few groups have carried out all of the proposed runs, most have done at least the baseline set of 
experiments defined in Table 1. NASA’s Global Modeling and Assimilation Office (GMAO) 
contributed runs made with version one of the NASA Seasonal-to-Interannual Prediction Project 
(NSIPP-1) AGCM. NOAA’s Climate Prediction Center, with support from the Climate Test 
Bed, contributed runs made with the Global Forecast System (GFS) AGCM, and NOAA’s 
Geophysical Fluid Dynamics Laboratory (GFDL) contributed runs made with the AM2. 1 
AGCM. The Lamont-Doherty Earth Observatory of Columbia University contributed runs made 
with the NCAR CCM3.0 AGCM, and NCAR contributed runs made with the CAM3.5 AGCM. 
An additional set of runs was made by COLA/University of Miami with the coupled 
(atmosphere-ocean) CCSM3.0 model employing a novel adjustment technique to nudge the 
coupled model towards the imposed SST forcing patterns. The main characteristics of the 
models of interest and some of the relevant references are presented in Table 2. 


Table 2. Description of the models used in the drought working group simulations. 


Model 

Resolution 

Convection Scheme 

Land Surface Model 

AM2pl2 

Delwoith at al. 2006 
GFDL Development Team 
(2004) 

2° x 2.5°, L24 

Relaxed Arakawa-Schubert 
(Moorthi and Suarez 1992) 

Milly and Shmakin 
(2002) 

GFS v2 

Campana, K. and P. 
Caplau (2005) 

T62 (~2°x2°), L64 

Simplified Arakawa 
-Schubert (Grell 1993; Pan 
and Wu 1995) 

Ek et al. (2003) 

mu 

3° x 3.75°, L34 

Relaxed Arakawa-Schubert 
(Moorthi and Suarez 1992) 

Mosaic 

(Koster and Suarez 1996) 

CCM3.0 

Kiehl et al. (1998) 
Seager et al. (2005) 

T42 (~2.8° x 2.8°) with 18 
hybrid sigma levels 

Zhang and McFarlane (1995) 
Hack (1994) 

Bonan (1996) 

CAM3.5 

http://www.ccsm.ucar.edu 

/models/atm-cam/ 

T85 with 27 hybrid sigma 
levels 

Oleson et al. (2008) 

Community Land Model 
Oleson et al. (2008) 
Stockli et al. (2007) 
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CCSM3.0 

AGCM: T85 with 26 levels 

Zhang and McFarlane (1995) 

Community Land Model 

Collins et al. 2006 

OGCM: l°xl° telescoping to 




l°xl/2° in deep tropics with 40 




levels 




3. Results 

This section provides ail overview of the results from the five AGCMs. Results from the 
coupled model (CCSM3.0) will be reported on in a separate paper. The main focus here is on 
the annual mean response over the United States to the Pacific SST anomaly pattern. Additional 
diagnostics are presented that summarize the results from all 8 combinations of the Pacific and 
Atlantic forcings shown in Table 1, as well as provide some assessment of the seasonality of the 
responses. The section ends with a brief overview of the responses to the trend pattern. 

We begin by examining the ability of the AGCMs to reproduce the observed annual mean 
precipitation and height field climatologies based on the available AMIP-style simulations (runs 
forced with observed SSTs) from each model for the period 1980-1998. We note that the results 
from the AMIP runs are quite similar to those from the control runs ( PnAn - see Table 1, not 
shown). The AMIP runs, however, provide a cleaner more direct comparison with observations 
than the control run, and therefore provide a more useful baseline assessment of model 
performance. 

All the models produce quite reasonable annual mean stationary wave and precipitation patterns 
(Figure 2). There are, however, differences in the details including such features as the strength 
of the Pacific ITCZ (c.f. the NSIPP1 and GFS results), and of particular interest here, the 
strength and orientation of the height anomalies over the Pacific/North American region. In 
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general, the models tend to overestimate precipitation in the western Pacific and Atlantic warm 
pool regions compared with observations. Most models also tend to underestimate the strength 
and eastward extent of the east Asian trough. The GFS model is the exception showing a deeper 
trough together with a stronger ridge over North America compared with the observations and 
the other models. An interesting feature is the wave train that appears to emanate from the 
central tropical Pacific and extends across North America into the Atlantic. It appears to be 
responsible for modifying the structure and amplitude of the west coast ridge and the trough over 
eastern North America. This feature is evident in the observations and all the models except 
GFS, where it is at most very weak. The extent to which such differences in the climatological 
stationary wave pattern impact the response to the SST anomalies is unclear. This issue will be 
revisited later in the discussions of the model responses. 

Figure 3 provides an overview of the agreement between the models with respect to the annual 
mean responses to the eight combinations of the Pacific and Atlantic SST patterns. The results 
are presented as spatial correlations between the 10 different combinations of the 5 models 
(model 1 correlated with model 2, model 1 correlated with model 3, etc.). The scatter of the 10 
different combinations gives a sense of the full range of agreement or disagreement between the 
various models for any one forcing pattern. The ordinate is the correlation based on 
precipitation, and the abscissa is the correlation based on the 200mb height. Figure 3a shows the 
results for the global and annual mean distributions. Here we limit the comparisons to the 4 
basic individual Pacific and Atlantic forcing patterns (no combinations of forcing patterns) to 
better highlight the differences in the level of agreement between responses to the Pacific and 
Atlantic forcing patterns. The global annual mean results show a number of interesting features. 
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First, it is clear that the models are in much stronger agreement regarding the responses to the 
Pacific forcing (dark colors), compared with the response to the Atlantic forcing (light colors). It 
is also clear that in general there is more agreement in the height response than in the 
precipitation response (almost all points lie below the 45° line). For plots 3b-f the results are 
shown for all 8 forcing patterns. The tropical results (Fig 3b) show generally very strong 
agreement in the height response for all forcing combinations. The response to the Atlantic 
forcing shows relatively weak precipitation correlations (generally less than 0.6), while the 
response to the combined (Pacific and Atlantic) forcing patterns exhibit correlations that are 
comparable to the Pacific forcing results - a reflection of the dominance of the Pacific forcing. 
There are also some interesting seasonal differences in the global correlations. The results for 
DJF (Fig. 3c) show a tendency for greater agreement in the precipitation response than the height 
response for the cold Pacific , while the opposite is true for the warm Pacific cases. This seems 
to reflect a weaker and less robust height response to a cold Pacific compared with a warm 
Pacific forcing. The height response appears to be most robust during MAM (Fig. 3d), while 
JJA shows overall the smallest precipitation correlations. The response to the Atlantic appears to 
be the most robust (and comparable to the response to the Pacific) during SON (3f). 

Figure 4 provides a more in-depth comparison of the annual mean global 200mb height and 
precipitation responses to the warm Pacific (PwAn). The results show considerable large-scale 
similarity among the models (as expected from the correlations in Figure 3). In particular, all 
models show the well-known (ENSO-type) horseshoe-shaped precipitation response with a 
positive precipitation anomaly in the central tropical Pacific surrounded by negative anomalies 
on either side of the equator. All the models also show negative anomalies over central 


13 



America, northeastern South America, and the tropical Atlantic. There is general agreement in 
the height anomalies with wave trains emanating from the tropical Pacific and extending 
poleward into both hemispheres. The precipitation anomalies over North America appear 
strongly coupled to the detailed structure of the height anomalies. In particular, the spatial extent 
of the positive precipitation anomalies over the United States is linked to the orientation and 
strength of the negative height anomalies over the continent. The response over the U.S. will be 
discussed further in reference to Figure 6. 

Figure 5 shows the global distribution of the responses to the cold Pacific (Pc An). The 
precipitation responses are largely of opposite sign compared with the responses to the PwAn 
pattern. In particular, all models again show the familiar horseshoe-shaped precipitation 
response with now a negative precipitation anomaly in the central tropical Pacific surrounded by 
positive anomalies on either side of the equator. All the models also show positive anomalies 
over Central America, northeastern South America, and the tropical Atlantic. The height 
anomalies also tend to be of opposite sign, but in the extratropics they tend to be weaker than the 
response to the warm Pacific, reflecting an asymmetric response to the warm and cold Pacific 
forcings. All the model responses show a tendency for a split ITCZ (also evident in the PwAn 
response) that is most pronounced in the GFDL and CAM3.5 model. All the models also show a 
ridge over the North Pacific that extends eastward across the United States (though less so for 
the GFS model). 
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Figure 6 provides a close-up of the U.S. surface temperature 4 and precipitation responses to the 
PwAn and PcAn SST patterns. While all the models show a tendency for wet conditions in 
response to PwAn forcing and dry conditions for the PcAn forcing, there are considerable 
differences among the models. For example, the NSIPP1 response tends to be relatively 
localized over the central Great Plains. On the other hand, the largest GFS response occurs 
along the southern and western tier of states, while the other models tend to show more 
widespread precipitation anomalies. The surface temperature and precipitation anomalies tend to 
correspond in the sense that wet anomalies are associated with cold anomalies and dry anomalies 
are associated with warm anomalies. The temperature response of the GFS model is quite 
different from the other models, showing extensive wanning over much of the US (cooling is 
confined to the southwest) in response to the PwAn forcing, and only a slight wanning over the 
southern Great Plains in response to the PcAn forcing. 

The linkage between the surface temperaUire and precipitation responses over the central United 
States is explored further in Figure 7. All the models show a very robust precipitation response 
with a clear separation of the wet (PwAn) and dry (PcAn) years. There is an overall tendency 
for a negative relationship between the precipitation and surface temperature anomalies. The 
exception is again the GFS model, for, which the response to the PwAn pattern is generally 
characterized by positive temperature anomalies and enhanced precipitation. The negative 
relationship between the annual mean temperature and precipitation responses is largely a 
reflection of the warm season responses (primarily JJA but also in some cases MAM and SON, 
not shown), and presumably reflects a strong tie between the atmosphere and land surface during 

4 Here and elsewhere in the text the surface temperature over land refers to the skin temperature. 


15 



the warm season for all but the GFS model 5 . Simply put, wetter conditions tend to lead to higher 
evaporative cooling and thus cooler air temperatures. 

Figure 8 summarizes the annual mean responses over the continental U.S. to the 8 different 
combinations of the Pacific and Atlantic forcing patterns. Figure 8a shows the precipitation 
responses. There is agreement among the models that a cold Pacific results in a reduction of 
precipitation, whereas a warm Pacific produces precipitation increases. Also, there is general 
agreement that a warm Atlantic leads to reduced precipitation whereas a cold Atlantic leads to 
increased precipitation, though with substantially smaller amplitudes, especially for the response 
to the cold Atlantic. Overall, the models agree that the combination of a cold Pacific and warm 
Atlantic (Pc Aw) tends to produces the largest precipitation deficits, whereas the combination of a 
warm Pacific and cold Atlantic (PwAc) tends to produce the largest precipitation surpluses. 

There is somewhat less agreement for the surface temperature responses (Figure 8b), with 
generally positive (negative) temperature anomalies associated with precipitation deficits 
(surpluses), but that is not the case for the GFS model (and to a lesser degree the NSIPP1 
model), which shows strong warm anomalies associated with enhanced precipitation (consistent 
with the previous results). 

The scatter in the year-to-year responses to the Pacific forcing shown in Figure 7 for the Great 
Plains shows the extent to which internal weather variability obscures the signal forced by the 
SST anomalies. This can be formalized in terms of the quantity 
R = (x-y)'s xy , 

1 The GFS model is in fact known to have a relatively weak land-atmosphere coupling strength (Koster et al. 2006). 
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where X and y represent seasonal values from the experiment and control runs respectively, and 


the overbar denotes a 50 (35 for the GFS model) year mean. Also, S 


2 

xy 



where 


£ v and S are the variance estimates of X and y, respectively. The mean difference between 

the experiment and control (the numerator) is a measure of the signal, so that R measures the 
size of the signal in units of standard deviation and can be viewed as the signal to noise ratio 
associated with the response. 


Figure 10 shows R for the precipitation response to the PwAn and PcAn forcing patterns for 
each season and all models over different regions of the U.S. (see Figure 9 for the definitions of 
the regions). The results for the U.S. average (top left panel) show significant (at the 5% level) 
responses throughout most of the year for both the PwAn and PcAn forcing patterns. Somewhat 
surprisingly the signal to noise ratio is smallest (in fact only marginally significant) in DJF and 
largest during MAM/AMJ, when it reaches values of 1 .5 or greater. We note that a very similar 
seasonality is found in the signal itself, so the results do not simply reflect a greater noise during 
the winter season. We caution however that, as indicated in Section. 2, there is no seasonal cycle 
in the prescribed forcing, whereas observed SST anomalies are largest in winter. 


A comparison of the results for PwAn and PcAn shows that magnitudes are largely comparable 
with perhaps somewhat weaker R values (characterized by less agreement among the models) 
for the cold Pacific case. The results for the U.S. as a whole reflect, to a large extent, the results 
in the Great Plains and SW, particularly the southern Great Plains. The northern Great Plains 
show a similar seasonality, but more modest signal to noise ratios. The NW U.S. has the largest 
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ratios during Northern spring and summer, while there is little evidence for a significant 
response in the NE. The SE has the largest signal to noise ratios during winter and spring with 
minimal (not significant) responses during the summer season. In general the models have a 
very similar seasonality of the response. The GFS model shows a somewhat different behavior 
especially for the northern Great Plains where the significant response is largely confined to the 
cold phase during the Fall. Also, in the SE the GFS model differs from the others in that it 
shows little significant response to the warm phase (PwAn). 

Figure 1 1 shows that the R values for the surface temperature response are substantially 
different from those for precipitation. In addition to having generally smaller amplitudes, there 
is considerable disagreement among the models regarding the seasonality and even the sign of 
R. For the continental U.S. as a whole (top left panel) the response to the warm Pacific is 
significant for 3 models during the warm season (May - August), and large and significant for 
one model (GFS) during the cold season (March - May). The latter result reflects the unusual 
(compared to the other models) southward extension of the upper level high in the GFS response 
(e.g.. Figure 4), and can be traced to the large signal-to-noise ratios in the northern tier of states 
(the second, third and fourth top panels of Figure 11). The response to the cold Pacific is for the 
most part not significant or marginally significant, with the largest values occurring during late 
summer and fall. Excluding the GFS model, the most consistent results and largest R values 
occur for the southern tier of states (bottom panels of Figure 1 1). For example, the southern 
Great Plains show a consistent warming signal during the warm season in response to a cold 
Pacific (with R peaking in late spring), and cold anomalies in response to a warm Pacific 
beginning in late winter and extending into late summer. The southwest shows considerable 
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asymmetry in the cold and warm Pacific responses, with only marginally significant responses to 
a cold Pacific whereas the response to the warm Pacific shows significant R values for two of 
the models during late spring and early summer. In the southeast (bottom right panel of Figure 
1 1) all the models agree on having little skill during the late summer to early winter seasons in 
response to the cold Pacific. Four models show a tendency for a significant warming in late 
spring and early summer in response to the cold Pacific , while three models show cooling during 
the summer in response to the warm Pacific. 

Finally, we present a brief overview of the responses of the models to the trend pattern. The 
focus is on the surface temperature response to the positive trend pattern and the results are 
compared to those from the companion set of runs in which the models were forced with a 
globally uniform SST wanning of 0. 16°C. The precipitation response to the trend (not shown) is 
weak for all the models, with anomalies that exceed 0.5 mm/day in amplitude largely confined to 
a few locations in the tropical Pacific and Indian Oceans. The basic response of all the models 
(left panels of Figure 12) is a tendency for warming over most of the world’s land areas, though 
there are substantial regional variations that differ between models. The GFS and GFDL models 
show the strongest warming, with especially large values (exceeding 0.5 °C) concentrated over 
North America and parts of Asia and Australia. CAM3.5 and CCM3 show the weakest 
warming, while the NSIPP1 model shows intermediate values with the largest warming 
occurring over North America. Focusing on North America, the GFS model shows strong 
wanning over most of the continent (especially the western half), while the GFDL and NSIPP1 
models show wanning that is more confined to the central and eastern United States. All but the 
GFS model show a substantial area of cooling spanning much of northern Canada and Alaska. 
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The impact of any regional variations in the SST Trend forcing pattern can be deduced from 
comparisons with the runs forced by the globally uniform wanning pattern 6 (right panels of 
Figure 12). The comparison shows that many of the features of the response to the trend pattern 
are reproduced in the response to uniform wanning. For example, some of the enhanced 
wanning over North America in the GFS model, the enhanced warming over Asia in the GFDL 
model, and more generally, the spatial pattern of the warming over much of Asia in all the 
models, is reproduced in the uniform warming case. The global spatial correlations between the 
responses to the trend and uniform wanning of the annual mean values over land range from 
0.34 for the NSIPP model to 0.43 for the GFDL model. There is, however, a large seasonal 
variation in the correlations for some models. For example, for the GFDL model, the 
correlations ranging from 0.7 in March to 0.15 in July. A key difference between the two sets of 
responses over North America is that the cooling in Northern Canada and Alaska noted earlier 
does not show up in the response to the uniform warming, indicating that this feature is primarily 
the result of regional variations in the SST trend pattern. Also, the localized warming responses 
over the U.S. in the GFDL and NSIPP models are not reproduced in the case of uniform 
warming. 

4. Summary and Discussion 

The USCLIVAR drought working group recently initiated a coordinated (multi-institutional and 
multi-model) effort to produce a set of idealized simulations designed to address fundamental 

6 This run was not done with CAM3.5. 
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questions regarding the physical mechanisms that link SST variations to regional drought, including 
an assessment of the role of land- atmosphere feedbacks. The set of experiments consist of 
multiyear simulations in which the models were forced by a number of idealized SST forcing 
patterns consisting of the leading rotated EOFs of SST variability on interannual and longer time 
scales. The main set of EOF forcing patterns include a global trend pattern, a Pacific ENSO-like 
pattern, and an Atlantic pattern that resembles the Atlantic Multi-Decadal Oscillation. Additional 
SST forcing patterns were designed to isolate ENSO and longer (decadal) time scales, and to isolate 
the influence of the tropical SST. A number of groups also ran experiments in which the land- 
atmosphere interactions were disabled by prescribing the soil moisture. 

This paper, in addition to providing a general overview of the project, attempts to provide a 
broad-ranged assessment of the model results focusing on overall behavior and highlighting 
where the models tend to agree and disagree. Results are limited to the responses to the two 
leading Pacific and Atlantic forcing patterns. While showing some aspects of the global-scale 
response, much of the focus of this paper is on the responses over the continental United States. 

A number of key results emerge from this initial analysis of the experiments. First, all the 
models produce similar (though different in detail) precipitation anomalies over the continental 
United States in response to the Pacific forcing pattern, with a tendency for reduced precipitation 
when forced with a cold Pacific and a tendency for enhanced precipitation when forced with a 
warm Pacific. The response to the Atlantic pattern is not as robust as the response to the Pacific , 
though there is a tendency for reduced precipitation when forced with a -warm Atlantic and a 
tendency for enhanced precipitation when forced with a cold Atlantic. There is general 
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agreement among the models that the largest precipitation response over the continental United 
States tends to occur when the two oceans have anomalies of opposite sign. That is, a cold 
Pacific and warm Atlantic tend to produce the largest precipitation reductions, whereas a warm 
Pacific and cold Atlantic tend to produce the greatest precipitation enhancements. 

The models tend to agree less on the area mean U.S. surface temperature response to both the 
Atlantic and Pacific forcing, though (with the exception of the GFS model) there is a general 
tendency for wet conditions to be associated with cold surface temperature anomalies and dry 
conditions to be associated with warm surface temperature anomalies. This greater disagreement 
for the surface temperature response appears to be the result of the sensitivity of the response 
over the U.S. to small (on planetary scales) shifts in the upper level wave response to the SST 
forcing. This is less so for the precipitation response. It is not clear exactly why that is the case, 
though it may be partly that the area mean of the more localized precipitation anomalies is less 
sensitive to shifts in the planetary wave forcing. The model differences are highlighted by the 
GFS response to the warm Pacific SST forcing, which tends to place the upper level ridge of the 
response over North America considerably further south compared to the other models, leading 
to warm surface temperature anomalies that extend well into the southern tier of states. 

Another key area of agreement among the models is in the seasonality of the signal-to-noise ratio 
( R ) of continental U.S. precipitation associated with the Pacific forcing. All the models show 
that the largest R values occur in spring with surprisingly small (not significant) R values 
during winter. The above results for the U.S. as a whole, reflect those in the Great Plains and the 
southwest, particularly the southern Great Plains. In contrast to these results, the R values of 
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the surface temperature response to the Pacific forcing are generally lower and show 
considerably less agreement among the models. For the continental U.S. as a whole the GFS 
model stands out as having very high R values during the cold season (March - May). This 
again appears to reflect the unusual (compared to the other models) southward extension of the 
upper level high in the GFS response and can be traced to the large R values that occur in the 
northern tier of states. For the other models, the most consistent results and largest R values 
occur for the southern tier of states, with, for example, the southern Great Plains showing a 
consistent warm anomaly during the warm season in response to a cold Pacific (with R peaking 
in late spring), and cold anomalies in response to a warm Pacific beginning in late winter and 
extending into late summer. 

The surface temperature response to the positive SST global trend forcing pattern shows 
substantial regional variations that are in part reproduced in runs forced with a globally uniform 
SST trend forcing. There is however substantial disagreement among the models in the 
regionality (e.g., the enhanced surface temperature response produced over North America by 
some of the models), highlighting the challenge of predicting regional impacts of global 
warming. The precipitation response to the trend forcing was found to be weak in all the 
models. 

The differences in the responses to the Pacific forcing pattern over the U.S. between the GFS 
and the other models is intriguing. Understanding these differences is important in view of the 
critical role that the GFS model (as part of the NOAA Climate Forecast System) plays in 
seasonal prediction, and more generally the increasingly important role that climate models play 
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in providing information on the regional impacts of global climate variability and change. While 
it is beyond the scope of this paper to fully address the reasons for these differences, it appears 
likely that the differences in the stationary wave patterns play a role. The GFS model produces a 
very reasonable stationary wave pattern and arguably produces the most realistic overall 
structure of the Pacific trough and North American ridge (Figure 2). The North American ridge 
is, however, stronger and the trough to the east considerably weaker compared with the other 
models and the observations. A preliminary linear model analysis of the NSIPP1 model 
stationary waves (results not shown) indicates that the eastern North American trough is 
particularly sensitive to the heating in the western tropical Pacific, suggesting that the 
differences we see in the climatological precipitation in that region may play an important role. 
On the other hand, most of the other models produce a rather weak Pacific trough - a problem 
that very likely also contributes to deficiencies in the response over North America. Clearly 
more work is needed to not only improve our understanding of the sensitivity of stationary 
waves to the climatological forcing, but also to determine the extent to which that sensitivity 
translates into uncertainties in the extratropical response to SST anomalies. 

The results of this initial analysis of the model experiments serve to highlight and quantify the 
important role of SST anomalies (especially those in the Pacific) in generating drought and 
pluvial conditions over the United States. The differences between the model results provide an 
assessment of the current uncertainties in our ability to model the global response to SST forcing 
(including the feedbacks associated with land-atmosphere interactions) and reinforce the need to 
improve our climate models. We expect that these results, as well as those reported in the other 
contributions to this special issue on drought will serve to stimulate further analysis of the 
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simulations, as well as suggest new research on the physical mechanisms contributing to 
hydroclimatic variability and change throughout the world. 
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List of Figures 


Figure 1: The three leading rotated EOFS (left panels) and associated PCs (right panels) of the 
annual mean SST based on the period 1901-2004. The values are scaled so that the product 
of the PCs and EOFs gives units of °C. The Pacific and Atlantic patterns reflect the two 
standard deviation forcing amplitude applied to the models, while the trend pattern must be 
divided by a factor of two to obtain the forcing amplitude. 

Figure 2: Annual mean precipitation and 200mb eddy height field averaged over the years 
1980-1998. The model results are from AMIP-style runs from each model (runs forced by 
observed SSTs). For the observations (lower right panel), the precipitation is from GPCP 
(Adler et al. 2003) and the height fields are from the NCEP/NCAR reanalyses (Kalnay et al. 
1996). The time period was chosen in order to have a common set of years for each model. 
Contour interval for the height field is 20m (negative values are dashed and the zero line is 
the first solid contour). Precipitation is in mm/day 

Figure 3: Spatial correlations of precipitation anomalies (ordinate) versus spatial correlations of 
200mb height anomalies (abscissa) for the various combinations of models and for the 8 
combinations of forcing for the Pacific and Atlantic patterns, a) global correlations of annual 
means, b) tropical correlations (+/- 30° latitude) of annual means, c) global correlations for 
DJF d) global correlations for MAM, e) global correlations for JJA, and f) global correlations 
for SON. See text for details. The colors refer to the forcing patterns indicated along the 
bottom of the plots. 
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Figure 4: Annual mean 200mb height and precipitation responses to the of PwAn SST anomaly 
pattern. The anomalies for each model are computed with respect to its annual mean response 
to a simulation in which the model is forced with climatological SSTs. See text for details on 
the SST forcing. Contour interval for height is 10 meters (negative values are dashed and the 
zero line is the first solid contour). Unit for precipitation is nim/day. 

Figure 5: Annual mean 200mb height and precipitation responses to the of PcAn SST anomaly 
pattern. The anomalies for each model are computed with respect to its annual mean 
response to a simulation in which the model is forced with climatological SSTs. See text 
for details on the SST forcing. Contour interval for height is 10 meters (negative values are 
dashed and the zero line is the first solid contour). Unit for precipitation is mm/day. 

Figure 6: Annual mean surface temperature (top panels, units: °C/day) and precipitation 

(bottom panels, units: mm/day) responses to the PwAn (left panels) and PcAn (right panels) 
SST anomaly pattern. The anomalies for each model are computed with respect to its 
annual mean response to a simulation in which the model is forced with climatological 
SSTs. See text for details on the SST forcing. 

Figure 7: Scatter plots of the Great Plains annual mean precipitation (ordinate) versus surface 
temperature (abscissa) response to the PwAn (red dots) and PcAn (blue dots) SST anomaly 
patterns for the Great Plains area average. Each point represents one of 50 (35 for the GFS 
model) years of each run. 


35 



Figure 8: The annual and continental United States mean responses for a) precipitation and b) 
surface temperature for all 8 combinations of the Pacific and Atlantic patterns for the 5 
AGCMs (see Table 1). 

Figure 9: The regions of the United States used to form the averages in Figures 10 and 11. 

Figure 10: Seasonality of the signal to noise ratios (R -see text) of the 3-month mean 
precipitation responses for each model over various regions of the United States to the 
Pacific warm (red curves) and cold (blue curves) SST anomaly pattern. The numbers along 
the abscissa refer to the center month of the 3-month means. Results are based on 50-year 
simulations except for the GFS model (dashed lines), which was run for 35 years. See text 
for the definition of the signal- to-noise ratio. The thin horizontal lines denote the 5% 
significance levels based on a t-test (for the GFS model the critical t-value is 0.49). Units: 
dimensionless. 

Figure 11: Seasonality of the signal to noise ratios (R -see text) of the 3-month mean surface 
temperature responses for each model over various regions of the United States to the 
Pacific warm (red curves) and cold (blue curves) SST anomaly pattern. The numbers along 
the abscissa refer to the center month of the 3-month means. Results are based on 50-year 
simulations except for the GFS model (dashed lines), which was mil for 35 years. See text 
for the definition of the signal- to-noise ratio. The thin horizontal lines denote the 5% 
significance levels based on a t-test (for the GFS model the critical t-value is 0.49). Units: 
dimensionless. 
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Figure 12: Left panels: The annual mean surface temperature response to the positive phase of 
the trend pattern. Right panels: The annual mean surface temperature response to a globally 
uniform SST wanning of 0.16 °C. The uniform wanning run was not performed with 
CAM3.5. Units: °C. 
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Figure 1: The three leading rotated EOFS (left panels) and associated PCs (right panels) 
of the annual mean SST based on the period 1901-2004. The values are scaled so 
that the product of the PCs and EOFs gives units of °C. The Pacific and Atlantic 
patterns reflect the two standard deviation forcing amplitude applied to the models, 
while the trend pattern must be divided by a factor of two to obtain the forcing 
amplitude. 
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Figure 2: Annual mean precipitation and 200mb eddy height field averaged over the years 1980-1998. 
The model results are from AMIP-style runs from each model (runs forced by observed SSTs). For the 
observations (lower right panel), the precipitation is from GPCP (Adler et al. 2003) and the height fields 
are from the NCEP/NCAR reanalyses (Kalnay et al. 1996). The time period was chosen in order to 
have a common set of years for each model. Contour interval for the height field is 20m (negative 
values are dashed and the zero line is the first solid contour). Precipitation is in mm/day. 
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Figure 3: Spatial correlations of precipitation anomalies (ordinate) versus spatial correlations of 
200mb height anomalies (abscissa) for the various combinations of models and for the various 
combinations of the Pacific and Atlantic forcing patterns, a) global correlations of annual means 
limited here to the 4 individual Pacific and Atlantic forcing patterns, b) tropical correlations (+/- 30° 
latitude) of annual means, c) global correlations for DJF d) global correlations for MAM, e) global 
correlations for JJA, and f) global correlations for SON. See text for details. The colors refer to the 
forcing patterns indicated along the bottom of the plots. 
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Figure 4: Annual mean 200mb height and precipitation responses to the of PwAn SST anomaly pattern. 
The anomalies for each model are computed with respect to its annual mean response to a simulation in 
which the model is forced with climatological SSTs. See text for details on the SST forcing. Contour 
interval for height is 10 meters (negative values are dashed and the zero line is the first solid contour). 
Unit for precipitation is mm/day. 
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Annual Mean precip_z200 (ColdPac-Clim) 
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Figure 5: Annual mean 200mb height and precipitation responses to the of Pc An SST anomaly pattern. 
The anomalies for each model are computed with respect to its annual mean response to a simulation in 
which the model is forced with climatological SSTs. See text for details on the SST forcing. Contour 
interval for height is 10 meters (negative values are dashed and the zero line is the first solid contour). 
Unit for precipitation is mm/day. 
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Figure 6: Annual mean surface temperature (top panels, units: °C/day) and precipitation (bottom panels, 
units: nnn/day) responses to the PwAn (left panels) and Pc An (right panels) SST anomaly pattern. The 
anomalies for each model are computed with respect to its annual mean response to a simulation in 
which the model is forced with climatological SSTs. See text for details on the SST forcing. 
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Figure 7: Scatter plots of the Great Plains annual mean precipitation (ordinate) versus surface 
temperature (abscissa) response to the PwAn (red dots) and PcAn (blue dots) SST anomaly patterns for 
the Great Plains area average. Each point represents one of 50 (35 for the GFS model) years of each 
run. 
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Figure 8: The annual and continental United States mean responses for a) precipitation and b) surface 
temperature for all 8 combinations of the Pacific and Atlantic patterns for the 5 AGCMs (see Table 1). 
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U.S. Sub-regions 



Figure 9: The regions of the United States used to form the averages in Figures 10 and 11. 
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Seasonality of Signal/Noise for Precipjand Response over US subregions: Pac-Clim 

3-Mon RunWean 
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Figure 10: Seasonality of the signal to noise ratios (R -see text) of the 3-month mean precipitation 
responses for each model over various regions of the United States to the Pacific warm (red curves) and 
cold (blue curves) SST anomaly pattern. The numbers along the abscissa refer to the center month of 
the 3-month means. Results are based on 30-year simulations except for the GFS model (dashed lines), 
which was run for 33 years. See text for the definition of the signal-to-noise ratio. The thin horizontal 
lines denote the 3% significance levels based on a t-test (for the GFS model the critical t-value is 0.49). 
Units: dimensionless. 
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Seasonality of Signal/Noise for SfcTempJand Response over US subregions; Poc-Clim 
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Figure 11: Seasonality of the signal to noise ratios (R -see text) of the 3-month mean surface 
temperature responses for each model over various regions of the United States to the Pacific warm (red 
curves) and cold (blue curves) SST anomaly pattern. The numbers along the abscissa refer to the center 
month of the 3-month means. Results are based on 50-year simulations except for the GFS model 
(dashed lines), which was run for 35 years. See text for the definition of the signal-to-noise ratio. The 
thin horizontal lines denote the 5% significance levels based on a t-test (for the GFS model the critical 
t- value is 0.49). Units: dimensionless. 
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Figure 12: Left panels: The annual mean surface temperature response to the positive phase of the trend 
pattern. Right panels: The annual mean surface temperature response to a globally uniform SST warming 
of 0.16 °C. The uniform warming run was not performed with CAM3.5. Units: °C. 
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Appendix A: Auxiliary Experiments 


In addition to forcing the models with the three main SST patterns presented in the text (the Pacific , 
Atlantic, and trend), the participating groups were encouraged to force their models with other 
patterns consisting of the tropical-only version of the Pacific and Atlantic patterns, and low and high 
frequency versions of the Pacific SST patterns. Details about these patterns are discussed below. 

In order to differentiate between the impacts of the tropics and extra tropics, another set of SST 
“tropical” forcing patterns was produced based on the interannual Pacific and Atlantic patterns. For 
the Pacific case, this was done by linearly tapering the values to zero between 15° and 21° latitude. 
The taper is such that the full amplitude occurs at 15°, Vi of the full amplitude remains at 18°, and the 
anomaly is identically zero from 21° latitude to the pole. Also, the meridional edges are such that 
the western boundary occurs at 120°E. For the Atlantic case, the Atlantic pattern was modified so 
that the edges of the box with the full anomalies were chosen as 88° W to 13°W, and 12°N to 18°N. 
The anomalies were tapered linearly north and south, with latitudes 9°N and 21°N getting 1/2 the 
anomaly, and with the anomaly going to 0 at latitudes 6°N and 24°N. 

In addition to the patterns described above (based on annual mean SST), two other patterns were 
produced for the Pacific that attempt to separate the ENSO and longer-term patterns of variability. 
The long time scales were isolated by applying a filter to the monthly SST data that retains time 
scales of about 6 years and longer (Zhang et al. 1997). The high (residual) frequencies (shorter than 
6 years) were obtained by subtracting the low frequency filtered data from the unfiltered monthly 
data. The leading REOFs and associated PCs from both the low pass (time scales of 6 years and 
longer) and residual (time scales shorter than 6 years) data are shown in Figure Al. In the case of 
the low-pass we focus on the second REOF (the first is again the trend pattern shown in Figure 1). 
The second low frequency REOF shows the well-known meridionally extensive pan-Pacific decadal 
pattern of variability (e.g., Barlow et al. 2001), with substantial middle latitude amplitude that in the 
Northern Hemisphere is 
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Figure Al: Top panels: The second REOF and PC of the low pass filtered (time scales greater than 6 
years) monthly SST data. Bottom panels: The first REOF and PC of the high pass filtered (time 
scales less than 6 years) monthly SST data. The results are based on the period 1901-2004. The 
values are scaled so that the product of the REOF and PC gives units of °C. 


linked to the Pacific Decadal Oscillation (Zhang et al 1997). In contrast, the leading REOF of the residual 
(high frequency) SST shows a clear ENSO structure with significant amplitudes that are largely confined to 
the central and eastern equatorial Pacific. 
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A final set of experiments was designed to assess the impact of soil moisture feedbacks. In this case 
archives of standard multi-decadal AMIP-style simulations were used to derive the climatological seasonal 
cycles of soil moisture content for each soil layer at each land point on the globe. Here weekly resolution 
was preferred, but in cases where these were not available monthly data were used. Subsets of the above 
idealized SST simulations were then repeated in such a way that at each time step, the simulated soil 
moisture states were thrown out and replaced with states interpolated from the climatological values 
established in the first simulation. 

Details of the various experiments and other information and links relevant to the 
USCLIVAR drought working group project may be found at: 

http://gmao.gsfc.nasa.gov/research/clivar drought wg/index.html 
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